GRASS GIS 101

This notebook is a quick introduction to GRASS and Jupyter Notebooks.

import sys
v = sys.version_info
print(f"We are using Python {v.major}.{v.minor}.{v.micro}")

We are using Python 3.10.12

GRASS Project Setup

# import standard Python packages
import os
import subprocess
import json
from io import StringIO
from pathlib import Path
from IPython.display import IFrame

Setup GRASS environment

sys.path.append(
    subprocess.check_output(["grass", "--config", "python_path"], text=True).strip()
)

# import GRASS GIS python packages
%reload_ext autoreload
%autoreload 2
import grass.script as gs
import grass.jupyter as gj

# create a temporary folder where to place our GRASS project
import tempfile
tempdir = tempfile.TemporaryDirectory()

Create a new GRASS project

gs.create_project(path=tempdir.name, name="foss4g2024", epsg="2817", overwrite=True)

# start GRASS in the recently created project
session = gj.init(Path(tempdir.name,"foss4g2024"))

# !ogrinfo "/vsizip//vsicurl/https://www.stlouis-mo.gov/data/upload/data-files/stl_boundary.zip" -al
# !ogrinfo /vsizip//vsicurl/https://www.stlouis-mo.gov/data/upload/data-files/streets.zip -al

Import data

City Boundary

gs.run_command("v.import",
   input="/vsizip//vsicurl/https://www.stlouis-mo.gov/data/upload/data-files/stl_boundary.zip",
   output="stl_boundary",
   layer="stl_boundary",
   snap=0.0001,
   # epsg=8801
)

Streets

streets_url="/vsizip//vsicurl/https://www.stlouis-mo.gov/data/upload/data-files/streets.zip"

gs.run_command("v.import",
   input=streets_url,
   output="stl_streets",
   layer="TgrGeoCd",
   snap=0.0001,
   # epsg=8801
)

map = gj.InteractiveMap(map_backend="folium")
map.add_vector("stl_boundary", fill_color="none", line_color="red")
# map.add_vector("stl_streets", line_color="grey")
map.add_layer_control()
map.show()

gs.run_command("g.region", vector="stl_boundary", res=10, flags="ap")

projection: 99 (NAD83(HARN) / Missouri West)
zone:       0
datum:      nad83harn
ellipsoid:  grs80
north:      298260
south:      271030
west:       1214160
east:       1226670
nsres:      10
ewres:      10
rows:       2723
cols:       1251
cells:      3406473

https://stacindex.org/catalogs

gs.run_command("t.stac.catalog", url="https://earth-search.aws.element84.com/v1/", format="plain", flags="")

---------------------------------------------------------------------------
Catalog: Earth Search by Element 84
---------------------------------------------------------------------------
Client Id: earth-search-aws
Client Description: A STAC API of public datasets on AWS
Client STAC Extensions: []
Client catalog_type: ABSOLUTE_PUBLISHED
---------------------------------------------------------------------------
Collections: 9
---------------------------------------------------------------------------
Collection: Sentinel-2 Pre-Collection 1 Level-2A 
---------------------------------------------------------------------------
Collection Id: sentinel-2-pre-c1-l2a
Sentinel-2 Pre-Collection 1 Level-2A (baseline < 05.00), with data and metadata matching collection sentinel-2-c1-l2a
Extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['2015-06-27T10:25:31.456000Z', None]]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Collection: Copernicus DEM GLO-30
---------------------------------------------------------------------------
Collection Id: cop-dem-glo-30
The Copernicus DEM is a Digital Surface Model (DSM) which represents the surface of the Earth including buildings, infrastructure and vegetation. GLO-30 Public provides limited worldwide coverage at 30 meters because a small subset of tiles covering specific countries are not yet released to the public by the Copernicus Programme.
Extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['2021-04-22T00:00:00Z', '2021-04-22T00:00:00Z']]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Collection: NAIP: National Agriculture Imagery Program
---------------------------------------------------------------------------
Collection Id: naip
The [National Agriculture Imagery Program](https://www.fsa.usda.gov/programs-and-services/aerial-photography/imagery-programs/naip-imagery/) (NAIP) provides U.S.-wide, high-resolution aerial imagery, with four spectral bands (R, G, B, IR).  NAIP is administered by the [Aerial Field Photography Office](https://www.fsa.usda.gov/programs-and-services/aerial-photography/) (AFPO) within the [US Department of Agriculture](https://www.usda.gov/) (USDA).  Data are captured at least once every three years for each state.  This dataset represents NAIP data from 2010-present, in [cloud-optimized GeoTIFF](https://www.cogeo.org/) format.

Extent: {'spatial': {'bbox': [[-160, 17, -67, 50]]}, 'temporal': {'interval': [['2010-01-01T00:00:00Z', '2022-12-31T00:00:00Z']]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Collection: Copernicus DEM GLO-90
---------------------------------------------------------------------------
Collection Id: cop-dem-glo-90
The Copernicus DEM is a Digital Surface Model (DSM) which represents the surface of the Earth including buildings, infrastructure and vegetation. GLO-90 provides worldwide coverage at 90 meters.
Extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['2021-04-22T00:00:00Z', '2021-04-22T00:00:00Z']]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Collection: Landsat Collection 2 Level-2
---------------------------------------------------------------------------
Collection Id: landsat-c2-l2
Atmospherically corrected global Landsat Collection 2 Level-2 data from the Thematic Mapper (TM) onboard Landsat 4 and 5, the Enhanced Thematic Mapper Plus (ETM+) onboard Landsat 7, and the Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) onboard Landsat 8 and 9.
Extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['1982-08-22T00:00:00Z', None]]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Collection: Sentinel-2 Level-2A
---------------------------------------------------------------------------
Collection Id: sentinel-2-l2a
Global Sentinel-2 data from the Multispectral Instrument (MSI) onboard Sentinel-2
Extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['2015-06-27T10:25:31.456000Z', None]]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Collection: Sentinel-2 Level-1C
---------------------------------------------------------------------------
Collection Id: sentinel-2-l1c
Global Sentinel-2 data from the Multispectral Instrument (MSI) onboard Sentinel-2
Extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['2015-06-27T10:25:31.456000Z', None]]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Collection: Sentinel-2 Collection 1 Level-2A
---------------------------------------------------------------------------
Collection Id: sentinel-2-c1-l2a
Sentinel-2 Collection 1 Level-2A, data from the Multispectral Instrument (MSI) onboard Sentinel-2
Extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['2015-06-27T10:25:31.456000Z', None]]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Collection: Sentinel-1 Level-1C Ground Range Detected (GRD)
---------------------------------------------------------------------------
Collection Id: sentinel-1-grd
Sentinel-1 is a pair of Synthetic Aperture Radar (SAR) imaging satellites launched in 2014 and 2016 by the European Space Agency (ESA). Their 6 day revisit cycle and ability to observe through clouds makes this dataset perfect for sea and land monitoring, emergency response due to environmental disasters, and economic applications. This dataset represents the global Sentinel-1 GRD archive, from beginning to the present, converted to cloud-optimized GeoTIFF format.
Extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['2014-10-10T00:28:21Z', None]]}}
License: proprietary
---------------------------------------------------------------------------
---------------------------------------------------------------------------
Conforms To:
---------------------------------------------------------------------------
    https://api.stacspec.org/v1.0.0/core
    https://api.stacspec.org/v1.0.0/collections
    https://api.stacspec.org/v1.0.0/ogcapi-features
    https://api.stacspec.org/v1.0.0/item-search
    https://api.stacspec.org/v1.0.0/ogcapi-features#fields
    https://api.stacspec.org/v1.0.0/ogcapi-features#sort
    https://api.stacspec.org/v1.0.0/ogcapi-features#query
    https://api.stacspec.org/v1.0.0/item-search#fields
    https://api.stacspec.org/v1.0.0/item-search#sort
    https://api.stacspec.org/v1.0.0/item-search#query
    https://api.stacspec.org/v0.3.0/aggregation
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/core
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/oas30
    http://www.opengis.net/spec/ogcapi-features-1/1.0/conf/geojson
---------------------------------------------------------------------------
---------------------------------------------------------------------------


gs.run_command("t.stac.collection",
               url="https://earth-search.aws.element84.com/v1/",
               collection="cop-dem-glo-30",
               format="plain",
               flags=""
            )

#     type: Collection
#     id: cop-dem-glo-30
#     stac_version: 1.0.0
#     description: The Copernicus DEM is a Digital Surface Model (DSM) which represents the surface of the Earth including buildings, infrastructure and vegetation. GLO-30 Public provides limited worldwide coverage at 30 meters because a small subset of tiles covering specific countries are not yet released to the public by the Copernicus Programme.
---------------------------------------------------------------------------
    links:
#         rel: self
#         href: https://earth-search.aws.element84.com/v1/collections/cop-dem-glo-30
#         type: application/json
#         rel: license
#         href: https://spacedata.copernicus.eu/documents/20126/0/CSCDA_ESA_Mission-specific+Annex.pdf
#         title: Copernicus DEM License
#         rel: parent
#         href: https://earth-search.aws.element84.com/v1
#         type: application/json
#         rel: root
#         href: https://earth-search.aws.element84.com/v1/
#         type: application/json
#         title: Earth Search by Element 84
#         rel: items
#         href: https://earth-search.aws.element84.com/v1/collections/cop-dem-glo-30/items
#         type: application/geo+json
#         rel: http://www.opengis.net/def/rel/ogc/1.0/queryables
#         href: https://earth-search.aws.element84.com/v1/collections/cop-dem-glo-30/queryables
#         type: application/schema+json
#         rel: aggregate
#         href: https://earth-search.aws.element84.com/v1/collections/cop-dem-glo-30/aggregate
#         type: application/json
#         method: GET
#         rel: aggregations
#         href: https://earth-search.aws.element84.com/v1/collections/cop-dem-glo-30/aggregations
#         type: application/json
---------------------------------------------------------------------------
    stac_extensions:
---------------------------------------------------------------------------

     item_assets:
---------------------------------------------------------------------------

         data:
#             type: image/tiff; application=geotiff; profile=cloud-optimized
---------------------------------------------------------------------------
            roles:
#             title: Data
---------------------------------------------------------------------------
            raster:bands:
#                 sampling: point
#                 data_type: float32
#                 spatial_resolution: 30
#                 unit: meter
#         data: {'type': 'image/tiff; application=geotiff; profile=cloud-optimized', 'roles': ['data'], 'title': 'Data', 'raster:bands': [{'sampling': 'point', 'data_type': 'float32', 'spatial_resolution': 30, 'unit': 'meter'}]}
#     item_assets: {'data': {'type': 'image/tiff; application=geotiff; profile=cloud-optimized', 'roles': ['data'], 'title': 'Data', 'raster:bands': [{'sampling': 'point', 'data_type': 'float32', 'spatial_resolution': 30, 'unit': 'meter'}]}}
#     title: Copernicus DEM GLO-30
---------------------------------------------------------------------------

     extent:
---------------------------------------------------------------------------

         spatial:
---------------------------------------------------------------------------
            bbox:
#         spatial: {'bbox': [[-180, -90, 180, 90]]}
---------------------------------------------------------------------------

         temporal:
---------------------------------------------------------------------------
            interval:
#         temporal: {'interval': [['2021-04-22T00:00:00Z', '2021-04-22T00:00:00Z']]}
#     extent: {'spatial': {'bbox': [[-180, -90, 180, 90]]}, 'temporal': {'interval': [['2021-04-22T00:00:00Z', '2021-04-22T00:00:00Z']]}}
#     license: proprietary
---------------------------------------------------------------------------
    keywords:
---------------------------------------------------------------------------
    providers:
#         name: European Space Agency
---------------------------------------------------------------------------
        roles:
#         url: https://spacedata.copernicus.eu/documents/20123/121286/CSCDA_ESA_Mission-specific+Annex_31_Oct_22.pdf
#         name: Sinergise
---------------------------------------------------------------------------
        roles:
#         url: https://registry.opendata.aws/copernicus-dem/
#         name: OpenTopography
---------------------------------------------------------------------------
        roles:
#         url: https://doi.org/10.5069/G9028PQB
#         name: AWS
---------------------------------------------------------------------------
        roles:
#         url: https://registry.opendata.aws/sentinel-1
#         name: Element 84
---------------------------------------------------------------------------
        roles:
#         url: https://element84.com
---------------------------------------------------------------------------

     summaries:
---------------------------------------------------------------------------
        gsd:
---------------------------------------------------------------------------
        platform:
---------------------------------------------------------------------------
        proj:epsg:
---------------------------------------------------------------------------
        storage:platform:
---------------------------------------------------------------------------
        storage:region:
---------------------------------------------------------------------------
        storage:requester_pays:
#     summaries: {'gsd': [30], 'platform': ['tandem-x'], 'proj:epsg': [4326], 'storage:platform': ['AWS'], 'storage:region': ['eu-central-1'], 'storage:requester_pays': [False]}

gs.run_command("t.stac.item",
               url="https://earth-search.aws.element84.com/v1/",
               collection="cop-dem-glo-30",
               format="plain",
               extent="region",
               method="bicubic",
               resolution="value",
               memory="5000",
               resolution_value=30,
               flags="d"
            )

Downloading Asset: {'href': 's3://copernicus-dem-30m/Copernicus_DSM_COG_10_N38_00_W091_00_DEM/Copernicus_DSM_COG_10_N38_00_W091_00_DEM.tif', 'type': 'image/tiff; application=geotiff; profile=cloud-optimized', 'title': 'Data', 'raster:bands': [{'sampling': 'point', 'data_type': 'float32', 'spatial_resolution': 30, 'unit': 'meter'}], 'roles': ['data'], 'collection_id': 'cop-dem-glo-30', 'item_id': 'Copernicus_DSM_COG_10_N38_00_W091_00_DEM', 'file_name': 'cop-dem-glo-30.Copernicus_DSM_COG_10_N38_00_W091_00_DEM.data', 'datetime': '2021-04-22T00:00:00Z'}
Import Url: /vsis3/copernicus-dem-30m/Copernicus_DSM_COG_10_N38_00_W091_00_DEM/Copernicus_DSM_COG_10_N38_00_W091_00_DEM.tif
Importing: cop-dem-glo-30.Copernicus_DSM_COG_10_N38_00_W091_00_DEM.data

gs.run_command("g.rename", raster="cop-dem-glo-30.Copernicus_DSM_COG_10_N38_00_W091_00_DEM.data,cop-dem-glo-30")

map = gj.InteractiveMap(width=500, height=500, map_backend="ipyleaflet", use_region=True)
map.add_raster("cop-dem-glo-30")
map.add_vector("stl_boundary", fill_color="none", line_color="red")
map.add_layer_control()
map.show()

map = gj.Map(width=300)
map.d_rast(map="cop-dem-glo-30")
map.d_vect(map="stl_boundary", color="red", fill_color="none")
map.show()

gs.run_command("g.region", raster="cop-dem-glo-30", res=30, flags="ap")
gs.run_command("r.colors", map="cop-dem-glo-30", color="elevation")
gs.run_command("r.slope.aspect", elevation="cop-dem-glo-30", slope="slope", aspect="aspect", overwrite=True)
gs.run_command("r.colors", map="slope", color="slope")
gs.run_command("r.relief", input="cop-dem-glo-30", output="relief", overwrite=True)

projection: 99 (NAD83(HARN) / Missouri West)
zone:       0
datum:      nad83harn
ellipsoid:  grs80
north:      298260
south:      271020
west:       1214160
east:       1226670
nsres:      30
ewres:      30
rows:       908
cols:       417
cells:      378636

map = gj.Map(width=300)
map.d_shade(shade="relief", color="cop-dem-glo-30")
map.d_vect(map="stl_streets", color="grey", fill_color="none", type="line", width=0.5)
map.d_vect(map="stl_boundary", color="red", fill_color="none")
map.show()

gs.run_command("r.geomorphon", elevation="cop-dem-glo-30", search=15, flat=0.1, forms="geomorphon", overwrite=True)
gs.run_command("r.mask", vector="stl_boundary")
map = gj.Map(width=300)
map.d_shade(shade="relief", color="geomorphon")
map.d_vect(map="stl_boundary", color="red", fill_color="none")
map.show()

gs.run_command("r.mask", flags="r")

gs.run_command("v.to.rast", input="stl_streets", output="stl_streets", use="dir", type="line", overwrite=True) 
gs.run_command("r.slope.direction", elevation="cop-dem-glo-30", direction="stl_streets", aspect="aspect", overwrite=True)

import ipywidgets as widgets
widgets.IntSlider()

from ipyleaflet import Map, Marker, basemaps, basemap_to_tiles
m = Map(
  basemap=basemap_to_tiles(
    basemaps.NASAGIBS.ModisTerraTrueColorCR, "2017-04-08"
  ),
  center=(52.204793, 360.121558),
  zoom=4
)
m.add_layer(Marker(location=(52.204793, 360.121558)))
m